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Abstract 

Nonstationary Gabor frames, recently introduced in adaptive signal anal- 
ysis, represent a natural generalization of classical Gabor frames by allowing 
for adaptivity of windows and lattice in either time or frequency. Due to the 
lack of a complete lattice structure, perfect reconstruction is in general not 
feasible from coefficients obtained from nonstationary Gabor frames. In this 
paper it is shown that for nonstationary Gabor frames that are related to some 
known frames for which dual frames can be computed, good approximate re- 
construction can be achieved by resorting to approximately dual frames. In 
particular, we give constructive examples for so-called almost painless nonsta- 
tionary frames, that is, frames that are closely related to nonstationary frames 
with compactly supported windows. The theoretical results are illustrated by 
concrete computational and numerical examples. 

Keywords: adaptive representations, nonorthogonal expansions, irregular 
Gabor frames, reconstruction, approximately dual frame 



1 Introduction 

Adapted and adaptive signal representation have received increasing interest over 
the past few years. As opposed to classical approaches such as the short-time Fourier 
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transform (STFT) or wavelet transform, adaptive representations allow for a variation 
of parameters such as window width or sampling density over time, frequency or both. 
Changing parameters in the frequency domain leads, for example, to non-uniform 
filter banks while adapting window width and sampling density in time is reminis- 
cent of the approach suggested in the construction of nonuniform lapped transforms. 
Transforms featuring simultaneous adaptivity in time and frequency are notoriously 
difficult to construct and implement, cp. fT5\ [7J [T2]; however they have shown to be 
useful in some applications , cf. [H]. On the other hand, fast and efficient imple- 
mentations exist for representations with adaptivity in only time or frequency. One 
recent method to obtain this kind of representations is represented by nonstationary 
Gabor frames, first suggested in [TT] and further developed in [U [TSJ [TO]. All the 
known implementations rely on compactness of the used analysis window in either 
time of frequency. This assumption allows for usage of tools developed for pain- 
less non-orthogonal expansions [6]. While a priori very convenient, the restriction to 
using compactly supported windows in the domain for which one wishes a flexible 
representation can be undesirable. As an example, we mention the construction of 
nonuniform filter banks via nonstationary Gabor frames, in which case this restriction 
forbids finite impulse response (FIR) filters; the latter are, however, imperative for 
real-time processing applications. 

In the current contribution, we therefore go beyond the results presented in the 
references above and consider nonstationary Gabor frames with fast decay but un- 
bounded support. The existence of this kind of frames was shown in [8]. Here we are 
concerned with methods for approximate reconstruction for these adaptive systems. 

2 Notation and Preliminaries 

Given a non-zero function g G L 2 (M), a modulation, or frequency shift, operator Mm 
is defined by Mug{t) : = e 27Tlblt g(t), and time shift operator T a k by T a kg(t) := g{t — ak). 
A composition, g^ i = MuT a kg(t) := e 2nMt g(t — ak) is a time-frequency shift operator. 

The set G(g, a, b) = {gu.i '■ k, I 6 Z} is called a Gabor system for any real, positive 
a, b. Q(g, a, b) is a Gabor frame for L 2 (R), if there exist frame bounds < A < B < oo 
such that for every / e L 2 (IR) we have 

A\\f\\l<J2\(L9k,i)\ 2 <B\\f\\l (1) 

k,iez 
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When working with irregular grids, we assume that the sampling points form 
a separated set: a set of sampling points {a k : k G Z} is called 5-separated, if 
\a>k — a>m\ > S ior a k , a m , whenever k 7^ m. xi wm denote the characteristic function 
of the interval /. 

A convenient class of window functions for time- frequency analysis on L 2 (M) is 
the Wiener space. 

Definition 1. A function g G L°°(M) belongs to the Wiener space ^(L 00 ,^ 1 ) if 

\\g\\w(L^^) ■■= X^ ess sn PteQ\9(t + k)\<oo, Q = [0, 1] . 

fcez 

For g G W(L°°, t 1 ) and 5 > we have [9j 

esssup t6K ^|#(t-5£;)| < (1 + (T 1 ) ||#|| W (L°°,£i) ■ (2) 

For / G L 2 (M) we use the following normalization of the Fourier transform, which 
we denote by T: 

JFf{u) = f(u) = [ f(t)e- 2 ^dt. 
Jr 

3 Nonstationary Gabor Frames 

Nonstationary Gabor systems are a generalization of classical Gabor systems of reg- 
ular time-frequency shifts of a single window function. 

Definition 2. Let g = {g k G W^L 00 ,^ 1 ) : k G Z} be a set of window functions and 
let b = {bk : k G Z} be a corresponding sequence of frequency-shift parameters. Set 
g k j = M hk ig k . Then, the set 

G(g t b) = {g k y. k,leZ} (3) 
is called a nonstationary Gabor (NSG) system. 

In generalization of regular Gabor frames, for which g k = T ak g, we will usually 
assume that the windows g k are localized around points a k in a separated set of time- 
sampling points {a k : k G Z}. Further, we will always make the assumption that the 
frequency sampling parameters b k are positive numbers contained in a closed interval, 
i.e. b k G [b L , b v ] C M + for all k G Z. 
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To every collection (J3J) we associate the analysis operator C g given by (C g f)k,i = 
if,9k,i), and synthesis operator U g , where U g c = J2k iez c k,Wk,i and c G £ 2 . For two 
Gabor systems £7(g, b) and £(7, b) the composition S 9il = U^C g , 

S ga f=J2(f^ M ^9k)M lbklk , (4) 

k,iez 

admits a Walnut representation for all / G £ 2 (M), [8]: 



S g «f(t) = Yl h ~k l 9k(t - lb- k l ) lk (t)f(t - lb~ k l ) . (5) 

k,iez 

We will frequently use the following correlation functions of a pair of Gabor systems: 
G!' 7 (t) = J2b^\g k (t - lb^)\\ lk (t)\, for I G Z. (6) 

k&L 

Note that this definition is asymmetric with respect to g and 7. Using this notation, 
we obtain the following bounds for the frame operator (j3J): 

II II 2 < ess sup ^^Gf' 7 • ess sup G]' g . (7) 
lez lez 

By inspection of (131), we note that the summands corresponding to / 7^ may be seen 
as the off-diagonal entries of the frame operator. We thus isolate the diagonal part 

Gr(t) = J2 b ^(t)\\lk(t)\ (8) 

fcez 

and denote the off-diagonal entries as follows: 

J R 9i7 = esssup Y S & fc 1 bfc(--% 1 )ll7fc(-)l- ( 9 ) 

iez\{o} k&L 

Note that, if g = 7, then the diagonal part of the frame operator S 9t9 is equal to G g ' g . 
Using this notation, we obtain the following additional bound: 

(S 9 , 7 f, f) WfWf < esssupG^ + y/Rg~~R^g , (10) 

Bessel sequences are of particular importance in the theory of frames and Riesz 
bases, [3 [19]. In the regular Gabor case, where g k (t) = g(t — ak) for some a > 0, it is 
sufficient to assume g G W(L°°,£ ) to obtain a Bessel sequence. We next provide a 
generalization of this property to NSG frames. 
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Proposition 3.1. Let Q(g,b) be a NSG system. If g k G W(L°°, i 1 ) for all k G Z 

with sup fceZ 1 1 j^fe 1 1 tv .e 1 ) bounded, and J2kez\9k(t)\ < B almost everywhere for some 
B < oo, then the sequence g k) i is a Bessel sequence. 

Proof. Let / G L 2 (IR). Then by the assumption on the windows g k and estimate ((?]) 

EKWI 2 = (S g *f,f) < 11/112 -esssup^Gf 9 

k,iez zez 

■ ess sup E&fc 1 !^^-) | \9k(- ~ ^fc 1 )! 



2 " 



Hln • ess sup |(1 + OlMlw^ 00 ,* 1 ) 



|2 



< II/II2 " 5 " sup[(l + b k )\\g k \\ w{L ^ 

fcez 



□ 



Given a frame, it is well known that there exists at least one dual frame G(-y, b) 
such that 

/= X)</,7fc,t>0M> forall/GL 2 (R). (11) 

The canonical dual frame is given by 7^ = S 1-1 ^/. In the regular Gabor case, 
the dual frames are again Gabor frames, i.e., they consist of time- frequency shifted 
versions of one dual window. This is due to the fact, that the frame operator S com- 
mutes with time-frequency shifts, hence j k ,i — S~ 1 9k,i = S~ 1 MuT ak g = MbiT ak S~ 1 g = 
MbiTak'y. In general, we cannot expect, that the dual frame of a NSG frame is again 
a NSG frame. 

However, even in the case of regular Gabor frames, it is often difficult to calculate a 
dual frame explicitly. For that reason, alternative possibilities to obtain perfect or ap- 
proximate reconstruction have been proposed, [U[2]. The following lemma quantifies 
the reconstruction error using general pairs of Bessel sequences. 

Lemma 3.2. Let Q(g,b) and Q(-y,h) be two Bessel sequences. Then 



5*9,7 1| < 



i-yXW/fc + y/Rg^ ■ R~f,g ■ (12) 

^ 00 

fcez 

Proof. Starting from the Walnut representation of S 9t7 , we estimate using Cauchy- 
Schwartz inequality for sums and integrals and, since all summands have absolute 
value, Fubini's theorem to justify changing the order of summation and integral: 
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i&z\{o} fcez 



< 



fcez 



ke2 



Rg^ ■ R^^g || J || 



since 



E ^ 1 7*(-w--^ 1 )/(--^ 1 )./ 



(13) 



(14) 



zez\{o} 



< E ^/bfcWIW-ifefc 1 )!!/^-^ 1 )!!/^)!* 



fee 
zez\{o} 



< 



E b ? [ / w - "rtiwoii/c* - ^)i 2 * 



£6Z\{0} 



1/2 



w-^ibikWii/wr* 



1/ 



< 



l/WI 2 E ^kWIW*-^ 1 )!^ 



fcez 
/ez\{o} 



1/2 r 



l/(*)l 2 E K l \i k mg k (t-iK l )\c 



fcez 
iez\{o} 



□ 



A special class of NSG systems are collections of compactly supported windows. 
They were first addressed in pp. The collection Q(g,b) with windows g k being com- 
pactly supported with | supp ^ | < 5- for all k is a frame for L 2 (M) if there exist 
constants A > and B < oo such that 

A < G 9 ' 9 (t) < B a.e. . (15) 

In this situation, Q(g,h) is called painless NSG frame. The canonical dual atoms 
are given by 7^ = M^ fe (Go' s 'g)' 1 gk- Note again that, in general, we may have 7^ = 
S~ 1 (Mi bk g k ) 7^ Mi bk (S~ 1 gk). If = b for all fc, then the frame operator commutes 
with the frequency- shifts and the dual frame is an NSG frame. 

The existence of NSG frames with not necessarily compactly supported windows 
was established in [8]. For these frames, finding canonical dual frames requires the in- 
version of the frame operator. This computation is expensive since the operator has 
considerably less structure than the frame operator in the classical, regular Gabor 
frame case, for which fast algorithms now exist, [TTf [131 KB] - To circumvent the prob- 
lem, we suggest the use of windows other than canonical duals to obtain sufficiently 
good approximate reconstruction. 
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4 Approximately dual atoms 

The notion of approximately dual pairs was discussed in jl]. For NSG Bessel sequences 
we adapt their definition as follows. 

Definition 3. Two Bessel sequences G(g, b) and £(7, b) are said to be approximately 
dual frames if \\I — S gn \\ < 1 or \\I — S J;g \\ < 1. 

Note that the two conditions given in the definition are equivalent since 

||i - s g j = \\i- c g u4 = \\i- u;c;\\ = \\i- c 7 u g \\ = \\i- s,j. 

In Definition [3] it is implicitly stated that, if two Bessel sequences are approximately 
dual frames, then each of them is a frame. This result was proved in [I] for general 
frames and it will be useful to reformulate the conditions for the existence of NSG 
frame, given in [S], in the context of approximately dual frames. 

Proposition 4.1. Let Q(g, b) be a Bessel sequence with Bessel bound B and < 
^1 — Y2kez\9k{t)\ 2 < A 2 < 00 a.e. for some positive constants Ai,A 2 . 

i) The multiplication operator Gq' 9 is invertible a.e. and, for 7^ = (Gq' 9 ) -1 ^, 

\\i-s 9 j< / ';%., ; . (i6) 

essmfGo 

ii) If 

R g , g < essinf Gg' 9 , (17) 
then Q(g, b) and ^(7, b) are approximately dual frames for L 2 (IR). 

Hi) Assume, additionally, for some 5-separated set of time-sampling points {ak ■ 
k G Z} and constants < Pu,Cl,Cjj < 00 such that for pk E]2,pu] C WL, 
Ck G [Cl,Cu] we have 

\g k {t)\ <C k (l + \t-a k \)- pk forallkeZ. (18) 

Then there exists a sequence {6°}fc 6 z ; such that for all sequence b k < b k , k G Z, 
(TTTD holds. 



Remark 1. If ( !T7|) holds, ^(7, b) is called a single preconditioning dual system for 
£(g,b). 
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Proof. Since all frequency modulation parameters bk are taken from a closed interval 
in IR + , the invertibility of Gq 9 is straightforward and the windows ik are well de- 
fined. Moreover, since G(g, b) is a Bessel sequence, so is G(~f, b), with Bessel bound 
(essinf Gq 9 )~ 2 B. Substituting 7^ = (G 9 ) ,9 )~ 1 gk for 7^ in the proof of Lemma |3~2| the 
first term in ( TT3|) vanishes and we obtain (i): 

K/-^/)i<|((G ! rr 1 E T / b k 1 9k(-)9k(--ibf)f(--ibf)j) 

zez\{o} kez 

< (essinf G^)^ E E^I^IW-^ 1 )!^--^ 1 )!.!- 

zez\{o} fcez 

< (essinf Gg-")" 1 ^!!/!!!. (19) 

(ii) follows directly from Definitional Finally, (Hi) follows from [HI Theorem 3.4], 
where it is shown that the assumptions (1T81) on the windows gu guarantee the existence 
of a sequence b° k such that 

(S ga f, /> II/II2 2 > essinf G 9 > 9 - R g , g > . (20) 

□ 

Single preconditioning dual windows are a good choice for reconstruction, when- 
ever the frame operator is close to diagonal. This is the case, if the original windows 
gk decay fast and frequency sampling is fast. 

If the frame of interest is close to some other frame, which, ideally, is better under- 
stood, other approximate dual windows may be derived from this frame. The proto- 
typical situation is a NSG frame which is close to a painless NSG frame in the sense of 
a small perturbation. Approximately dual frames in the context of perturbation the- 
ory were recently studied in [1]. In the following proposition we give error estimates 
for the reconstruction with approximately dual frames in such a situation. This pro- 
vides different reconstruction methods apart from single preconditioning which was 
addressed in Proposition 14.11 

Proposition 4.2. Assume that £(g,b) is a Bessel sequence with bound B and that 
Q(h, b) is a NSG frame with lower and upper frame bound Ah and Bh, respectively. 
We set vpk = hk — gt and define the following windows: 

( a ) Ik 1 = $h h,hk,i (canonical dual of hkj) (21) 

(b) lli =Sl\9k,i (22) 
Then the following hold: 
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(i) 

\\I-S g „4 <^ 1/2 ||^||. (23) 

If ess sup ^2 leZ Gf^ < Ah, then Q(^ l ,h) and (?(g, b) are approximately dual 
frames. 

(ii) 

\\I-S gr/ 2\\ <A-\^ h + VB)\\C4. (24) 

If ess sup Y^iez ^ < (^b~+Vb) ~ ' then Q{~{ 2 ,h) and £(g,b) are approximately 
dual frames. 

Remark 2. The first statement of Proposition 14.21 is contained in |4]. 
According to [8], the assumption that ess sup Gf'^ < A/j can be satisfied if 
the functions ipk decay polynomially, i.e., (V'fcWI — Cfc(i + |t|)~ Pfc with appropriate 
constants Ck and decay rates pk > 1 . 

Proof. From ([7]) it follows that ||C^|| 2 < esssup J2iez Gf^ ■ The same estimate holds 
for H^ll 2 . 

Since Q{h,h) is a frame with canonical dual frame ^(7 1 ,b), an upper frame bound 
of Gij 1 , b) is given by A^ 1 . We thus obtain (i) as follows: 

\\I-S g „i\\ = WU^Ch-U^CgW < ||Z7 7 i||||C^|| <^ 1/2 ||C V ,||. (25) 

If esssup J2iez Gf ^ < A h , then \\I — S g ^i \\ < 1 and (/(7 1 , b) and G(g, b) are approx- 
imately dual frames as claimed. 

To show (ii), we note that U^2 = S^ h U g and thus 

11-^ ~~ ^,7 2 H —W^h,h^h,h — S hth S g> g\\ = WSfJ^UhCh — UgCg)\\ (26) 
<\\S^ h \\\\U h C h - UgC h + UgC h - UgCg\\ = A^ \\U h ~ U gC ^\\ (27) 

<Al l \\C4(^ET h + VB) (28) 
and (IMD follows. □ 



4.1 Perturbation of painless nonstationary Gabor frames 

If a NSG system can be derived as a perturbation of a painless NSG frame, the approx- 
imately dual windows given in Proposition 14.21 are particularly simple to compute. In 
this situation, the frame Q(h, b) is the painless frame {?(g°, b) and the frame operator 
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Sh,h = S g o g o is the multiplication operator Gq ,g . Moreover, in this particular case, 
the approximately dual frames, given by 7^ = (Gq ' 9 )~ 1 g k i an d ill = (^0 9 )~ 1 9k,i 
are NSG frames. This is a very important asset, since for NSG frames fast algorithms 
for analysis and reconstruction using FFT exist. 

In [8] we constructed a special class of NSG frames, arising from painless NSG 
frames, which we introduce next. 

Definition 4 (Almost painless NSG frames). Let (?(g, b) be a NSG system, assume 
that the windows g k are essentially bounded away from zero on the intervals I k = 
[a k - (2b k y l ,a k + (2b k )~ 1 ] and set g° k = g k Xi k - If £(g°,b) is a (painless) frame for 
L 2 (IR), then we call the system Q(g,b) an almost painless NSG system (or frame). 

For almost painless NSG systems, the estimates given in Proposition 14.21 can be 
written more explicitly. 

Corollary 4.3. Assume that Q(g, b) is an almost painless NSG system and let A = 
essinf Gq ' 9 to be the lower frame bound of the painless frame £?(g°,b) and g T k = 
gk — gkXh ■ Then the following hold: 

(i) fori l k = {Gl°> 9 °)-'gl, 

\\I ~~ Sg^l || < Aq 1 a/ Rgr, g o ■ Rg° , g r (29) 

(n) fo ri l = (G^°)-'g k 

\\I ~~ ^g,i 2 \\ — 1 I Rg r ,g° + Rg°,g r + ess SU P 9 ) 

V i& ) 

Proof. The estimates follow from Lemma 13.21 First, substituting 7^ = (Gq ' 9 )~ l g k 
for 7fc in (TT2l . the first term vanishes, since ^2 ke zK 1 '9k9k = &o ' § > an( I we obtain 

\\I — S gr/ i\\ < \/ R g>7 i • R^,g < (essinf Gq ' 9 ) _1 a/ Rg r ,g° ' Rg°,g r j 

since 

^ 7l =esssup £ E^W-^IK^)" 1 ^-)! ( 31 ) 

zez\{o} fcez 

^(essinfGf^r-esssup £ ^^|^(- - ^ 1 )||^(-)| 

;ez\{o} k&L 

= (essinf Gq ' 9 ) _1 R g r j9 o , 
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similarly for i? 7 i 



7 ,S" 



For (ii) we substitute 7^ = {G 9 Q ,9 )~ 1 gk for 7 in the proof of Lemma [3^21 g k + gl for 
g k and use the fact that g° k and g r k have disjoint supports. Then, since g k (- — lb k l )g k {-) 
is zero for / 7^ 0, using (fT4l) we obtain that 



< 



fcez 



+ (essinf(G 



9°,9°\-l 



9' ,3' 



The first term vanishes and by ([7]), H^^r || < esssup^ ZgZ G 



g' ,g' 



□ 



5 Examples 

We present two examples to illustrate our theory. In the first example we deal with 
almost painless frames using Gaussian windows. We consider three different approx- 
imately dual frames and check their performance in terms of reconstruction. 

In both examples, we consider a basic window and dilations by 2 and |, respec- 
tively. Since the dilation parameters take only three different values, there are three 
kinds of windows, with support size 1/2, 1 and 2, respectively. Note that, while 
theoretically possible, sudden changes in the shape and width of adjacent windows 
turn out to be undesirable for applications, hence we only allow for stepwise change 
in dilation parameters. 

Example 5.1. Let s k G { — 1, 0, 1} with \s k — s k -i\ G {0, 1} for all fceZ. We consider 
a sequence of windows g k that are translated and dilated versions of the Gaussian 
window g(t) = e~ n[at]2 : g k (t) = T ak ^g(b k t) = y/hg(b k (t - a k )), with cr = 2.5, 
b k = 2 Sfc , a = and for all k G Z 

a k+1 = a k + (2b k y 1 if s k = s k+ i , 
a k+ i = a k + (36 fc+ i) _1 if s k > s k+1 , 
a k+1 = a k + (3b k )~ l if s k < s k+1 . 

Here, bi = 1/2, b\j = 2 and the {a k : k G Z} are separated with minimum distance 
S = 1/4. We arrange the windows as follows: after each change of window size, no 
change is allowed in the next step; in other words, each window has at least one 
neighbor of the same size. 
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Let Ik = [afc — (26fc) _1 ,afc + (26fc)" 1 ] and define a new set of windows by = 
9k(t)Xh- Then {M^g^ : fc, / £ Z} is a painless nonstationary Gabor frame with 
lower frame bound Aq = 0.1718. 

The system Q(g, b) arises from the painless frame £(g°, b), and therefore we are 
interested in the approximate dual windows proposed in Corollary 14.31 We first con- 
sider 7^ = (Gq ,9 the canonical dual frame of (?(g°,b). According to ( l29l) . we 
need to calculate R g r tg o and R g °, g r in order to obtain an estimate of the reconstruction 
error. 

Claim 1: R g r >g o < 0.00827 + C(10~ 6 ). 
For fixed k £ Z, 

zez\{o} lez 

since g r k and g k have disjoint support. Then, due to 6,7 ^periodicity of ^fc \9k(^~ 
Ibfr 1 )], we have 



R gr , ga = esssup^^ 1/2 |^(-)| E^ V2 |^(- " ^)l 

fcez ;ez 

< esssup^6- 1/2 |^(-)| ■ esssup t6lfe ^6" 1/2 |^(t - lb?)\ (32) 



fcez zez 



In order to obtain more accurate estimates, we split Ik by setting l£ = [a^, + 
(26fc)" 1 ] and 1% = [ak — (2bk)~ 1 ,ak\ and estimate expression (1321) by its values at the 
end points of l£ or 1^ , respectively. We observe that, for t £ I~£: 

oo 

iez z=i 

oo oo 

< ^ e -H 2 + ^ e -^-i)/2] 2 < 0073g + o(10 -6) 
;=1 1=1 

By the symmetry of g£ with respect to a^, an analogous estimate holds for t E I k ■ 

^ —1/2 

Furthermore, due to the arrangement of the windows g^, ess sup z2 keZ b k Ifl'fcOl — 

1.1206 and Claim 1 follows. 

Claim 2: R g o^ < 0.0157 + C(10" 6 ). 

Observe that E; e z\{o} K ^l^i* ~ l K V )\ < ll fo fe V2 ^IU = 1. Then, with J+ and 
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Figure 1: An example of an arrangement of windows 





1 1 
: 1 


I 

< >r< , 




1 1 




















m i ( m 












I ! i 
\ ' ' 
\ ' ' 


/l 


l 

\ 




■ 1 / 


V 

1 


V v , 






1 


2 


3 4 





Figure 2: A zoom-in on the tails g r k in the arrangement of Fig. [TJ 



I m defined as before: 



i? s o >g r = esssup^6 fc 1/2 |^(-)| ^ & fc 1/2 bfe(- -^fe 1 )! < esssup^6 fc 1/2 |^(-; 



zez\{o} 



feez 



< 



sup { ess sup t6J + 6 fe 1/2 |^(t) | , ess sup te/ - b k 1/2 l^fc W I } 

m£Z *> ' 



sup { esssup, e/ + J2 e ~ W[(Tbkit ' ak)f Xi^t) , esssup te/m ]Te^ CT 



— — 



c- 



We bound C + and C~ by their maximal values on I+, respectively I~. The set 
{a*; : A; G Z} is 5— separated, hence \a^ — a m \ > \k — m\S. By the arrangement of the 
windows, there are at most two windows g r k which assume their maximum e 



-7r[<x/2] 2 



111 



the interval or in I m . Without loss of generality we assume that the two maximal 
values occur in /+, corresponding to the windows g r m _i and g r m+2 - Then, g r m+1) g r m 
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are zero in /+, cf. Figure |2] for an example situation. Therefore, 

(J+ < e -TrW b k(am+(2b m )- 1 -a k )] 2 _|_ 2 e -^b k (2b k )- 1 ] 2 _|_ e -7r[o-6 fc (a m -a fe 

fc>m+2 k<m—l 

-7r[crfefcfc(5] 2 



< 2 + 

\ k>0;kb t 



e 

fc>2 



Since we assumed that two windows g r k reached their maximum in /+, C < C + . As 
before, the bound from Claim 2 follows by numerical calculations. 

In summary, due to ( 1291) . Claim 1 and 2 and the lower frame bound Aq = 0.1718, 
the reconstruction error using the approximate duals 7^ is bounded by 

\\I -S 9j7 i || < 0.0663 + e>(10~ 6 ). (33) 

Another choice of approximate dual system are the windows 7^ = (Gq ,9 )~ l gk- 
In this setting we use fl30|) to derive 

III-^H < A^ 1 l Rgr^o + R g o t9 r +esssup^Gf ' §r J < 0.1402 + C(10~ 6 ), 

since, by previous calculations, 

esssup^Gf * = esssup^6- 1/2 |^(t)| ^6~ 1/2 |^(t - lbf)\ (34) 

zez fcez zez 

< esssup^6~ 1/2 |^(t)| ■ esssup^fc" 172 !^^-^ 1 )! 

< 0.0001158 + C(10~ 6 ). 

As a third choice of approximate dual windows we consider single preconditioning 
windows 7^ = {G 9 Q 9 )~ 1 gk introduced in Proposition 14.11 It can easily be seen that 

Rg,g — Rg°,g r Rg T ,g° ^-g r ,g r — ^-g°,g r ^-g r ,g° ^SS SUp ^ ^ G 9 ,9 , (35) 

and, by previous calculations, it follows 

R g ,g < 0.0241 + C(10" 6 ) < A < ess inf G 9 ' 9 . (36) 

Therefore, by Proposition 14. II (ii) . Q(g, b) and £(7, b) are approximately dual frames. 
Moreover 
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Remark 3. Observe that from each of the approximately dual frame estimates given 
in the above example, the frame property of (?(g, b) follows. 

Note that, from ( 1341) . the frame property of G(g,b) may also be derived by apply- 
ing results from perturbation theory cf. [3]. Indeed, since i eZ \(f, gk,i — 9li)\ 2 = 
E k ,iez\(f,9 r k ,i)\ 2 < esssup£^Gf'' 5 l/ll! < A ||/||I, it follows that 0(g,b) is a 
frame with a lower frame bound A = 0.1630. 

Example 5.2. In our second example, we turn to the situation mentioned in the 
introduction, namely, the construction of non-uniform filter banks with compactly 
supported, that is, FIR filters, via NSG frames. In this situation, the frame operator 
S does not have a Walnut-like structure as given in ((5)) on the time side. However, 5" 
may be considered on the frequency side by applying a Fourier transform. Then, we 
encounter the same structure as before and may exploit the developed techniques to 
deduce the frame property and to construct approximate dual frames for reconstruc- 
tion. The situation is schematically depicted in Figure |3j 

Window arrangement Window arrangement on Fourier transform side 




Time Frequency 

Figure 3: An example for the arrangement of dilated windows in Example 2 

It is obvious that, while, in the time domain, various windows gk with different 
pass-bands, are shifted to obtain the overall system, applying a Fourier transform 
yields the known situation: various windows g^ are modulated to create an NSG 
system. 

On the other hand, we are interested in using FIR filters, that is, the windows g^ 
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are all compactly supported, hence not bandlimited. Similar to the construction in 
Example 15.11 we can cut the windows gk to obtain a painless NSG reference frame. 
More precisely, we consider the family of windows gk and a vector of corresponding 
time-shift parameters a^. Then, we set g^j = T ak igk and are interested in the frame 
property of the set of functions {gk,i '■ k,l G Z}. Considering, due to the lack of 
structure on the time side as mentioned above, the corresponding frame operator on 
the frequency side corresponds to investigating the operator FSJ 7 *, which acts on 
the Fourier transform of a signal of interest. In other words, we are now dealing with 
the set of functions {F(gk,i) = M ak i§k : k, I e Z}. 

For the current example, we consider Hanning windows hj. and, as in the previous 
example, apply dilations by 2 _1 and 2, respectively, to obtain various time- and 
frequency resolutions. The time-shift parameters are chosen in parallel to the 
choice of the frequency-shift parameters bk in Example 15.11 

Given the explicit knowledge of the spectral properties of the Hanning windows, 
explicit error estimates can be derived in a similar manner as in the previous exam- 
ple. Here, we also numerically calculate the errors resulting from reconstruction by 
means of the three different proposed approximate dual systems. We use the same 
nomenclature as before, that is, 7^, 7^ and 7)? denote the canonical duals of the pain- 
less frame, the set of windows (Gq ' 9 ) _1 (?fc and the single preconditioning windows, 
respectively. Then 



1. 


\\I - 


Sg,-/ 1 | 


= 0.0210 


2. 


\\I- 




= 0.0407 


3. 


|| J — 


Sg,y 3 I 


= 0.0407 



As before, the canonical duals of the painless frame provide the best approximate 
reconstruction. 

The set of windows §k used in this example, together with their true, single precon- 
ditioning duals 7^ and the canonical duals of the corresponding painless frame, 7^, 
are depicted in Figure HI On the right plots, zoom-ins are shown to better compare 
the detailed behavior. 

A comparison between a basic window from the true dual frame and the precondi- 
tioning duals 7^ and the true dual and the approximate dual 7^ is depicted in Figures [5] 
and O respectively, in both the frequency domain (upper plot) and the time domain 
(lower plot). It should be noted that, while compactly supported in frequency, 7^ 
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still have better decay in the time domain than both the true dual windows and the 
other approximate dual windows. 
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Original (almost painless) frame 



Original (almost painless) frame 




200 400 
True dual frame 




200 400 600 
Single preconditioning dual frame (y 3 ) 




200 400 
Canonical duals of painless frame (y 1 ) 
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True dual frame 




400 420 440 460 480 500 
Single preconditioning dual frame (y 3 ) 




400 420 440 460 480 500 
Canonical duals of painless frame (y 1 ) 
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Figure 4: The original windows and the windows from various (approximate) dual 
frames. 
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Solid: true dual, dotted: y 3 








200 250 300 350 400 




Figure 5: Comparison between true dual and single preconditioning dual 7 s . 



Solid: true dual, dotted:-/ 1 




Figure 6: Comparison between true dual and painless approximate dual 7 . 
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